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ABSTRACT 

Among the blazars detected by the Fermi satelUte, we have selected the 23 blazars that in the 
three months of survey had an average 7-ray luminosity above 10'*® erg s^^. For 17 out of 
the 23 sources we found and analysed X-ray and optical-UV data taken by the Swift satellite. 
With these data, implemented by archival and not simultaneous data, we construct the spectral 
energy distributions, and interpreted them with a simple one-zone, leptonic, synchrotron and 
inverse Compton model. When possible, we also compare different high energy states of sin- 
gle sources, like 0528+134 and 3C 454.3, for which multiple good sets of multi-wavelength 
data are available. In our powerful blazars the high energy emission always dominates the 
electromagnetic output, and the relatively low level of the synchrotron radiation often does 
not hide the accretion disk emission. We can then constrain the black hole mass and the disk 
luminosity. Both are large (i.e. masses equal or greater than IO^Mq and disk luminosities 
above 10% of Eddington). By modelling the non-thermal continuum we derive the power 
that the jet carries in the form of bulk motion of particles and fields. On average, the jet power 
is found to be slightly larger than the disk luminosity, and proportional to the mass accretion 
rate. 

Key words: BL Lacertae objects: general — quasars: general — radiation mechanisms: non- 
thermal — gamma-rays: theory — X-rays: general 



1 INTRODUCTION 

We would like to attack in a more systematic way than done in the 
past the problem of the relation between the power carried by the 
jet of blazars in the form of bulk motion of particles and magnetic 
fields, and of the luminosity associated to accretion. This task is 
not easy, mainly because of the Doppler enhancement of the non- 
thermal continuum produced by jets, and because this continuum 
very often hides the radiation produced by the accretion disk. But 
starting from the early work of Rawlings & Saunders (1991), ev- 
idence has been accumulated that the jet power of extragalactic 
radio-loud sources is at least of the same order as the accretion 
one, and possibly a factor ~10 larger. 

On the larger scale (radio-lobe size) one can resort to mini- 
mum energy considerations and estimates of the lifetime of radio- 
lobes to find out the minimum jet power needed to sustain the 
radio-lobe emission (Rawlings & Saunders 1991). The main uncer- 
tainty associated to this argument is the unknown energy contained 
in the proton component. 

At smaller, but still very large, jet scales (on kpc to Mpc), 
one can use the recently discovered (by the Chandra satellite) X- 
ray emission from resolved knots, to model it and to infer the to- 
tal number of leptons needed to produce the observed radiation. 
Here the main uncertainties are related to the number of low en- 
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ergy leptons (of energy 7min?TieC^ present in the source, a number 
that is rarely well constrained. Furthermore, one also needs to as- 
sume how many protons are associated for each emitting lepton. 
Assuming a one to one correspondence and using observed data 
to limit 7min, Tavecchio et al. (2000); Celotti, Ghisellini & Chi- 
aberge (2001), Tavecchio et al. (2004), Sambruna et al. (2006) and 
Tavecchio et al. (2007) derived jet powers that are comparable with 
those inferred at the sub-pc scales, suggesting that the jet powers 
in these sources are conserved all along the jet, with the caveat that 
for large scale jets is difficult to constrain the bulk Lorentz factor 
r. However, values compatible with the inner jet are preferred (see 
Ghisellini & Celotti 2001). 

A technique developed quite recently makes use of the cavities 
or "bubbles" in the X-ray emitting intra-cluster medium of cluster 
of galaxies, and measures the energy required to inflate such bub- 
bles. Assuming that this energy is furnished by the jet, one can 
calculate the associated jet power (i.e. Churazov et al. 2002; Allen 
et al. 2006; Balmaverde, Balbi & Capetti 2008). 

At the VLBI scale (pc or tens of pc), one takes advantage of 
the resolving power of VLBI to measure the size of the synchrotron 
emitting region: by requiring that the Self Compton process does 
not overproduce the total observed X-ray flux, one derives a lower 
limit to the Doppler factor 5 [defined as 5 = l/[r(l — P cos 6^)], 
where 9^ is the viewing angle and V the bulk Lorentz factor], and 
a limit on the number of emitting particles and on the value of the 
magnetic field (Celotti & Fabian 1993). This in turn gives an esti- 
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mate of the jet power. The main uncertainties here are the number 
of low energy electrons (they emit unobservable self-absorbed syn- 
chrotron emission) and the fraction of the total X-ray flux produced 
by the radio VLBI knot. 

The advent of the Compton Gamma Ray Observatory and its 
high energy EGRET instrument on one side, and of ground based 
Cherenkov telescopes on the other, finally let it possible to discover 
at what frequencies blazar jets emit most of their power, i.e. above 
MeV energies. Variability of the high energy flux also tells us that 
the emission site cannot be too distant from the black hole, while 
the absence of 7-7 absorption (leading to pair production) tells us 
that the emission site cannot be too close to the black hole and 
its accretion disk (Ghisellini & Madau 1996; Ghisellini & Tavec- 
chio 2009; hereafter GT09). Bracketed by these two limits, one ob- 
tains a few hundreds of Schwarzschild radii as the preferred jet 
location where most of the dissipation occurs. Modelling the ob- 
served spectral energy distribution (SED) from the mm to the 7- 
ray band returns the particle number and the field strength needed 
to account for the observed data (see e.g. Celotti & Ghisellini 2008, 
hereafter CG08; Kataoka et al. 2008, Maraschi & Tavecchio 2003). 
Again, the main uncertainty in the game is 7min, but in some cases 
7min can be strongly constrained. This occurs when the considered 
blazar has good soft X-ray data characterised by a flat spectrum. 
In these cases the X-ray emission is likely to be due to inverse 
Compton scatterings between low energy electrons and seed pho- 
tons produced externally to the jet (external Compton, EC; e.g. 
Sikora, Begelman & Rees 1994). In fact, while the synchrotron 
Self-Compton (SSC) process uses relatively high energy electrons 
and internal synchrotron radiation to produce X-rays, in the EC 
case we do see X-rays produced by electrons of very low energies, 
strongly constraining 7min (see e.g. Tavecchio et al. 2007). 

As for the accretion disk radiation, it is usually hidden by 
the stronger optical-UV synchrotron emission. But in very pow- 
erful blazars we have good reasons to believe that the synchrotron 
component originating the radio to IR-optical radiation is not com- 
pletely hiding the thermal emission produced by the accretion disk. 
These reasons come partly from past observations of relatively 
high redshifts blazars (even if not detected by Fermi or EGRET, 
see e.g. Landt et al. 2008; Maraschi et al. 2008; Sambruna et al. 
2007; GT09) and partly from theoretical considerations about the 
expected scaling of the magnetic energy density with the black hole 
mass: we expect lower magnetic fields in jets associated to larger 
black hole masses (e.g. GT09), and thus a reduced importance of 
the synchrotron component of the SED. Furthermore, the presence, 
in these objects, of broad emission lines of relatively "normal" 
equivalent width (compared to radio-quiet objects) ensures that the 
thermal ionising continuum cannot be much below the observed 
optical emission. 

The above considerations guide us to select the most powerful 
blazars as very good candidates for measuring both the jet power 
and the accretion luminosity. In turn, since in these sources most 
of the luminosity is emitted in the hard X-rays and in the 7-ray 
band, it is natural to take advantage of the recently published list of 
blazars detected by the Large Area Telescope (LAT) on board the 
Fermi Gamma Ray Space Telescope (Fermi). It revealed more than 
one hundred blazars with a significance larger than lOcr in the first 
three months of operation (Abdo et al. 2009a, hereafter A09). Of 
these, 57 are classified as flat spectrum radio quasars (FSRQs), 42 
as BL Lac objects, while for 5 sources the classification is uncer- 
tain. Including 2 radio-galaxies the total number of extragalactic 
sources amounts to 106. Redshifts are known for all FSRQs, for 30 
BL Lacs, for 1 source of uncertain classification and for the two 
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Table 1. The 23 most powerful Fermi blazar.s in the A09 catalogue. In 
the last 3 columns we indicate the logarithm of the average 7-ray lu- 
minosity as observed by Fermi during the first 3 months of survey (cgs 
units), if the source was detected by EGRET, if there are Swift observa- 
tions, and if there is the public available 7-ray light-curve at the web page: 
|http : / / f ermi . gsf c ■ nasa ■ gov/ ssc/data/access/lat/| . 

radio-galaxies, for a total of 90 objects. Within this sources, there 
are 23 blazars exceeding an average 7-ray luminosity of 10** erg 
s^^ within the 3 months survey. These are the targets of our study, 
aimed to estimate the power associated to accretion and to the jet. 
We will also take advantage of the optical-UV and X-ray obser- 
vations made by the Swift satellite, to construct the spectral energy 
distributions (SED) of Fermi blazars, a crucial ingredient to reli- 
ably constrain the emission models leading to the jet and accretion 
power estimate. 

In this paper we use a cosmology with h = JIa ~ 0.7 and 
f^M = 0.3, and use the notation Q = QxW'^ in cgs units (except 
for the black hole masses, measured in solar mass units). 



2 THE SAMPLE 

Tab. [T] lists 23 blazars taken from the A09 catalogue of Fermi de- 
tected blazars. We selected them simply on the basis of their K- 
corrected (see e.g. Ghisellini, Maraschi & Tavecchio 2009) 7-ray 
luminosity being larger than 10** erg s^^. All sources but two 
are Flat Spectrum Radio Quasars (FSRQs). The exceptions are AO 
0235+164 and PKS 0426-380, classified as BL Lac objects. How- 
ever, these sources do have broad emission lines, (see e.g. Raiteri et 
al. 2007 for AO 0235+164 and Sbarufatti et al. 2005 for PKS 0426- 
380), visible in low emission states. We therefore believe that AO 
0235+164 and PKS 0426-380 are FSRQs whose line emission is 
often swamped by the enhanced non-thermal continuum. In Tab. 
[T]we mark the 8 blazars for which there are public available Fermi 
light curves. Although we use for all blazars the 3-months averaged 
7-ray flux and spectral index, for these 8 sources we can check if at 
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the time of the Swift observations the 7-ray integrated photon flux 
was much different than the average (the spectral index is available 
only for the 3-months average). For all sources but PKS 1502+106 
we find a good consistency. For PKS 1502+106 the Swift observa- 
tions were probably performed during the rapid decay phase after a 
major 7-ray flare. Therefore it is difficult to assess the exact 7-ray 
flux during the Swift observation, but in any case it should not be 
different than the used average flux by more than a factor 2. 



3 SWIFT OBSERVATIONS AND ANALYSIS 

For 17 of the 23 blazars in our sample there are Swift observations, 
with the majority of the objects observed during the 3 months of the 
Fermi survey. The data were analysed with the most recent software 
Swift_Rel3 . 2 released as part of the Heasoft v. 6.6.2. 
The calibration database is that updated to April 10, 2009. The XRT 
data were processed with the standard procedures (XRTP IPELINE 
V . . 1 2 . 2). We considered photon counting (PC) mode data with 
the standard 0-12 grade selection. Source events were extracted in a 
circular region of aperture ~ 47", and background was estimated in 
a same sized circular region far from the source. Ancillary response 
files were created through the xrtmkarf task. The channels with 
energies below 0.2 keV and above 10 keV were excluded from the 
fit and the spectra were rebinned in energy so to have at least 30 
counts per bin. Each spectrum was analysed through XSPEC with 
an absorbed power-law with a fixed Galactic column density from 
Kalberla et al. (2005). The computed errors represent the 90% con- 
fidence interval on the spectral parameters. Tab.|2]reports the log of 
the observations and the results of the fitting the X-ray data with a 
simple power law model. 

UVOT (Roming et al. 2005) source counts were extracted 
from a circular region 5"— sized centred on the source position, 
while the background was extracted from a larger circular nearby 
source-free region. Data were integrated with the uvotimsum 
task and then analysed by using the uvotsouroe task. The ob- 
served magnitudes have been dereddened according to the formu- 
lae by Cardelli et al. (1989) and converted into fluxes by using stan- 
dard formulae and zero points from Poole et al. (2008). Tab.[3]list 
the observed magnitudes in the 6 filters of UVOT, and the Galactic 
extinction appropriate for each source. 



4 THE MODEL 

To interpret the overall SED of our sources we use a relatively 
simple, one zone, homogeneous synchrotron and Inverse Compton 
model, described in detail in GT09. This model aims at accounting 
in a simple yet accurate way the several contributions to the radia- 
tion energy density produced externally to the jet, and their depen- 
dence upon the distance of the emitting blob to the black hole. Here 
we summarise the main characteristics of the model. The emitting 
region, of size rdiss and moving with a bulk Lorentz factor F, is 
located at a distance Rdiss from the black hole of mass M. The ac- 
cretion disk emits a bolometric luminosity Ld. The jet is assumed 
to be accelerating in its inner parts, with V oc R^^'^ {R is the dis- 
tance from the black hole), up to a value Fmax. In the acceleration 
region the jet is assumed to be parabolic (following, e.g. Vlahakis 
& Konigl 2004); beyond this point the jet is conical with a semi- 
aperture angle tp (assumed to be 0.1 for all sources). 

The energy particle distribution Nl'y) [cm~^] is calculated 
solving the continuity equation where particle injection, radiative 
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Figure 1. SED of PKS 0048-071, 0202-17 and PKS 0215+015, together 
with the fitting models, with parameters listed in Tab. |4] Fermi and Swift 
data ai'e indicated by dark grey symbols (red in the electronic version), while 
archival data (from NED) are in light grey. The short-dashed line is the 
emission from the IR toras, the accretion disk and its X-ray corona; the 
long-dashed line is the SSC contribution and the dot-dashed line is the EC 
emission. The solid light grey line (green in the electronic version) is the 
non theiTnal flux produced by the jet, the solid dark grey line (blue in the 
electronic version) is the sum of the non-thermal and thermal components. 



cooling and pair production (via the 7-7 



process) are taken 



into account. The injection function (5(7) [cm s ] is assumed 
to be a smoothly joining broken power-law, with a slope (3(7) oc 



and 7 below and above a break energy 7b: 



Q(7) 



(7/7b)" 



1 + (7/7b)- 



(1) 



© 0000 RAS, MNRAS 000, 000-000 



4 G. Ghisellini, F. Tavecchio and G. Ghirlanda 



source 


OBS date 
dd/mm/yyyy 


^Gal 

1020 cm-2 


r 


x2/dof 


^0.2-10,unabB 
10-12 cgS 


^2-10,unabs 
10-12 CgS 




Zo/Uu/ZUUj 


Oft 
J. ZD 


1.4oitU.l 


1 II LJ 


1 1 _i_n 1 


(\-\-C\ 
Z.UltU.Z 


uzz /— joy 


U//1 1/ZUUo 


AO 
Z.4Z 


l.ZZltU. 14 


ft ft/7 
o.o/ / 


1 ft-l-H 


l.ZZltU. 1 J 




uz/uy/zuuo 


/. / 


1.44lllU. iZ 


on/i 4 

ZU/ 14 




1 ft-un 1 


A'J/1'7 oil 


i j/ iU/ZUUo 


A O 
4.Z 


l.jJitU.4/ 


0/1/1 
Z.4/J 


1 nj_n o< 


n 7J_n 1 
U. I±\j.j 


U4ZO— JoU 


0*7/1 n/onno 
Z //iU/ZUUiS 


O 1 

Z. i 


1 oi_i_n o 1 

i.yjitu.zi 


1/3 


1 <j_n 


n 7J_n 
U. /ztU.Z 


U4j4-Zj4 


OA/ 1 n/onno 
Zo/ iU/ZUUo 


O Q 
Z.iS 


i.DltU.j 


1/i 


1 oj_n o 
i.ZitU.Z 


U.iSitU.Z 


UjZ!5+i J4 


no /I n/onnQ 


O/l o 
Z4.Z 




1 Q//1 
1.5/4 


A n_Ln ft 




UiSZU+JuU 


on /I o/onno 
ZU/lZ/ZUUo 


4. / 1 


1 /;<j_n 1 ft 


1 QQ/ft 


1 1/1 j_n 1 ft 
1.04iLU.l0 


n Q-Un o 


uyzu+44 


lo/Ul/zUUy 


1 /IT 

1.4/ 


1. /4±U.lj 


11 /I /I o 

13.4/12 


A o_i_n A 
4.Z±U.4 


o /i_i_n o 
z.OrtU.z 


14D4— jD4i 


U//U1/ZUU(5 


D.D4 


1 QftJ_n Oyl 


1 C/y1 


1 i_i_n 1 


n ft_i_n 1 

U.DltU.lZ 


1454-3542 


12/09/2008 


6.54 


1.84±0.31 


1.06/1 


3.7±0.3 


1.8±0.4 


1502+106 


08/08/2008 


2.29 


1.55±0.10 


18.4/18 


2.6±0.2 


1.6±0.2 


1520+319 


12/11/2008 


2.0(*) 


1.17±0.4 


43/44 


0.6 




1633+382 


22/02/2007 




1.56±0.25 


7.7/5 


1.58 




2023-07 


08/12/2008 


3.37 


1.78±0.3 


3.8/3 


1.5±0.3 


0.7±0.3 


2227-088 


28/04/2005 


4.21 


1.55±0.12 


29.7/22 


4.5±0.4 


2.7±0.1 


2251+158 


08/08/2008 


6.58 


1.7±0.06 


485/527 


33±1 


518±1 


2325+093 


03/06/2008 


4.1 


1.11±0.17 


6/5 


3.7±0.6 


2.9± 0.5 



Table 2. Results of the X-ray analysis. *: poorly determined spectrum, the C-Statistic was used to fit the spectrum, a. the spectrum was fitted from 0.5 keV. 
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Table 3. Results of the UVOT analysis. The given magnitudes are not corrected for Galactic extinction. The value of Ay is taken from Schlegel et al. (1998) 



The total power injected into the source in the form of rel- 
ativistic electrons is P/ = me(?V f (3(7)7(^7, where V = 
(47r/3)rdiss is the volume of the emitting region, assumed spheri- 
cal. 

We assume that the injection process lasts for a light crossing 
time rdiss/c, and calculate N{'y) at this time. The rationale for this 
choice is that even if injection could last longer, adiabatic losses 
caused by the expansion of the source (that is travelling while emit- 
ting) and the corresponding decrease of the magnetic field would 
make the observed flux to decrease. Therefore our calculated spec- 
tra correspond to the maximum of a flaring episode. 

Besides its own synchrotron radiation, the jet can scatter pho- 
tons produced externally to the jet by the accretion disk itself (e.g. 
Dermer & Schlickeiser 1993), the broad line region (BLR; e.g. 



Sikora, Begelman & Rees 1994) and a dusty torus (see Blazejowski 
et al. 2000; Sikora et al. 2002). We assume that the BLR (assumed 
for simplicity to be a thin spherical shell) is located at a distance 
^?BLR = lO^^L^^g cm, emitting, in broad lines, a fraction /blr = 
0.1 of the disk luminosity. The square root dependence of Rblr 
with Ld implies that the radiation energy density of the broad line 
emission within the BLR is constant, but is seen amplified by a fac- 
tor ~ by the moving blob, as long as i?diss < Rblr- A dusty 
torus, located at a distance Rjr — 2.5 x lOi^L^^^ cm, reprocesses 
a fraction /ir (of the order of 0.1-0.3) of La through dust emission 
in the far IR. Above and below the accretion disk, in its inner parts, 
there is an X-ray emitting corona of luminosity Lx (we fix it at a 
level of 30% of La)- Its spectrum is a power-law of energy index 
ax = 1 ending with a exponential cut at Ec =150 keV. We cal- 
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Figure 8. SED of PKS 0528+134 at four epochs. For the first three, the 7-ray data comes from EGRET. For the 2008 SED, the 7-ray data are from Fermi. 
Symbols and lines as in Fig.[T] 



culate the specific energy density (i.e. as a function of frequency) 
of all these external components as seen in the comoving frame, to 
properly calculate the resulting Inverse Compton spectrum. 

5 FITTING GUIDELINES 

In this section we illustrate some considerations that are used as a 
guide for the choice of the parameters. 

• All our sources do have broad emission lines (including the 
two sources classified as BL Lacs, as discussed above). The pres- 
ence of visible broad lines ensures that the ionising continuum, that 
we associate to the accretion disk, cannot be much lower that the 
observed optical-UV flux. 

• The shape of the 7-ray emission, as indicated by the slope 
measured by Fermi, is related to the slope of the electron energy 
distribution in its high energy part. It is then associated to the slope 
of the high energy synchrotron emission. 

• The bolometric power of our sources is always dominated by 
the 7-ray flux. The ratio between the high to low energy emission 
humps (Lc/Ls) is directly related to the ratio between the radia- 
tion to magnetic energy density U^/U^. If the dissipation region is 



within the BLR, our assumption of -Rblr = 10^'' ^^45 cm gives 

where we have assumed that t/^' ~ t/BLR- 

• The peak of the high energy emission (i^c) is believed to be 
produced by the scattering of the line photons (mainly hydrogen 
Lyman-Q) with electrons at the break of the particle distribution 
(7pcak)- Its observed frequency, uc, is 

2 2 

VC ~ — I^LyaTS^pc^i, (3) 

I + z 

A steep (a > 1) spectrum indicates a peak at energies below 100 
MeV, and this constrains FSj^^^]^. For our sources, the scattering 
process leading to the formation of the high energy peak occurs 
always in the Thomson regime. 

• The strength of the SSC relative to the EC emission depends 
on the ratio between the synchrotron over the external radiation 
energy densities, as measured in the comoving frame, (7s'/^cxt- 
Within the BLR, U^xt depends only on P^, while Us depends on 
the injected power, the size of the emission, and the magnetic field. 
The larger the magnetic field, the larger the SSC component. The 
shape of the EC and SSC emission is different: besides the fact that 
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Figure 9. SED of PKS 2251-158=3C 454.3 at four epochs. For the July 2007 SED, 7-ray data comes from AGILE (Vercellone et al. 2009). For the 2008 SED, 
the 7-ray data are from Fermi. Symbols and lines as in Fig.[T] 



the seed photon distributions are different, we have that the flux at 
a given X-ray frequency is made by electron of very different ener- 
gies, thus belonging to a different part of the electron distribution. 
In this respect, the low frequency X-ray data of very hard X-ray 
spectra are the most constraining, since in these cases the (softer) 
SSC component must not exceed what observed. This limits the 
magnetic field, the injected power (as measured in the comoving 
frame) and the size. Conversely, a relatively soft spectrum (but still 
rising, in vF^) indicates a SSC origin, and this constrains the com- 
bination of B, rdiss and P/ even more. 

• In our powerful blazars the radiative cooling rate is almost 
complete, namely even low energy electrons cool in a dynamical 
timescale rdiss/c. We call 7cooi the random Lorentz factor of those 
electrons halving their energies in a timescale rdiss/c- Note that 
7cooi is rarely very close to unity, but more often a few). Therefore 
the corresponding emitting particle distribution is weakly depen- 
dent of the low energy spectral slope, si, of the injected electron 
distribution. 

• In several sources we see clear signs of thermal emission in the 
optical-UV part of the spectrum. If associated to the flux produced 
by a standard multi-colour accretion disk, we can estimate both the 
mass of the black hole and the total disk luminosity. The tempera- 



ture profile of a standard disk emitting locally as a black-body is 
(e.g. Frank, King &. Raine 2002) 



y4 _ 3RsLd 



1/2- 



(4) 



where Ld = riMc? is the bolometric disk luminosity and Rs is 
the Schwarzschild radius. The maximum temperature (and hence 
the peak of the uFi, disk spectrum) occurs ?lI R ^ 5Rs and 
scales as Tmax oc (Ld/iEdd)'^'^*Af The level of the emis- 
sion gives Ld [that of course scales as (Ld/iEdd) M]. Therefore, 
for a good optical-UV coverage, and when the synchrotron emis- 
sion is weaker than the thermal emission, we can derive both the 
black hole mass and the accretion luminosity. Given a value of the 
efficiency 77, this translates in the mass accretion rate. 

• The one-zone homogeneous model here adopted is aimed to 
explain the bulk of the emission, and necessarily requires a com- 
pact source, self-absorbed (for synchrotron) at ~ 10^^ Hz. The 
flux at radio frequencies must be produced further out in the jet. 
Radio data, therefore, are not directly constraining the model. In- 
directly though, they can suggest a sort of continuity between the 
level of the radio emission and what the model predicts at higher 
frequencies. 
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Table 4. List of parameters used to construct the theoretical SED. Subsciipts refer to the year of the corresponding SED or to an EGRET {E) state as opposed 
to a "Fermi" (F) state. Col. [1]: name; Col. [2]: redshift; Col. [3]: dissipation radius in units of 10^^ cm and (in parenthesis) in units of i?g; Col. [4]: black 
hole mass in solar masses; Col. [5]: size of the BLR in units of 10^^ cm; Col. [6]: power injected in the blob calculated in the comoving frame, in units of 10*^ 
erg s~^; Col. [7]: accretion dislc luminosity in units of 10^^ erg s~^ and (in parenthesis) in units of LEdd; Col. [8]: magnetic field in Gauss; Col. [9]: bulk 
Lorentz factor at Rdiss; Col. [10]: viewing angle in degrees; Col. [11] and [12]: break and maximum random Lorentz factors of the injected electrons; Col. 
[13] and [14]: slopes of the injected electron distribution [Q{-f)] below and above 7^,; For all cases the X-ray corona luminosity Lx = 0.3L^. Its spectral 
shape is assumed to be oc cxp{—hu/l50 kcV). 

6 RESULTS 

The SED of our sources, together with the best fitting model, are 
shown in Figs.[T]-[9] In Tab.|2]and Tab.[3]we list the result of our 
XRT and UVOT analysis for the 17 blazars with Swift data. In Tab. 
|4]we list the parameters used to compute the theoretical SEDs and 
in Tab.|5]we list the the power carried by the jet in the form of radi- 
ation, electrons, magnetic field and protons (assuming one proton 
per emitting electron). All optical-UV fluxes shown in Fig.[T|-|9] 
have been de-reddened according to the Galactic Av given in the 
NED database and reported in Tab. [3] but not for the (possible) ab- 
sorption in the host galaxy nor for Lya absorption (line and edge) 
due to intervening matter along the line of sight. According both 
to theoretical consideration (see IVIadau, Haardt & Rees 1999) and 
optical-UV spectra of high redshift quasars (e.g. Hook et al. 2003) 
this kind of absorption, in quasars, should be at most marginal. This 
is an important point, since it allows us to consider the bluer pho- 
tometric points of UVOT as a relatively good estimate of the in- 
trinsic flux of the source. Inspection of the optical-UV SED shown 
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in Figs.[T]-|9]strongly suggests that in some cases this emission is 
due to the accretion disk. This is when there is a peak in the i/Fi, 
optical-UV spectrum with an exponential decline. This cannot be 
the synchrotron peak, since the steep spectral index implied by the 
steep 7-ray energy emission (bow-ties) is incompatible with this 
interpretation. 



This fact, together with the theoretical considerations men- 
tioned in the previous section, allow us to consider the optical- 
UV emission of nearly half of our blazars as due to the accretion 
disk emission. Among the sources with UVOT data exceptions are 
0215+015, 0235+164, 0454-234, 0528+134, 1454-354, 2251-158 
and 2325+093. For these sources either the non-thermal continuum 
dominates the flux, or the quality of the data is not sufficient to dis- 
criminate between the thermal and the non-thermal processes. For 
them, the black hole mass and disk luminosity has been chosen 
(when summed to the non-thermal emission) not to overproduce 
the observed flux. 
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Log u [Hz] Log u [Hz] 

Figure 2. SED of PKS 0227-369; AO 0235+164 and PKS Q-iAl-lll. Sym- Figure 3. SED of PKS 0426-380, PKS 0445-234 and PKS 0820+560 Sym- 

bols and lines as in Fig.[T] bols and lines as in Fig.[T] 



6.1 Black hole masses, accretion and outflow rates 

The derived black hole masses are in the range (l-6)xlO^M0(Fig. 
Illb . and the disk luminosities, in Eddington units, cluster around 
Ld/ L-Edd ~ 0.1-0.3 (see Fig. llOl ). Fig.lUlshows the accretion and 
outflow mass rates (in units of Mq yr~^ and their ratio. All our 
sources requires more than 1 Mq yr~^ of accreting mass, while the 
outflow rate is a factor 0.01-0.1 smaller. This can be understood 
in simple terms, considering that in matter dominated jets Pj = 
FMoutC^ and then 

Maccr ^ i^d 1 1 Ld 



6.2 Jet power 

Table |5] lists the power carried by the jet in the form of radiation 
(Pr), magnetic field (Pb), electrons (Pc) and protons (Pp, assuming 
one proton per emitting electron). All the powers are calculated as 

P = nr^T^ficUi (6) 

where Ui is the energy density of the i component. 

The total power carried by the jet is often larger than the disk 
luminosity, and sometimes exceeds the Eddington luminosity (see 
Fig.llOland Fig.ll2t. In Fig. [12] we compare Pjct as a function of 
Ld for our sources with the FSRQs and BL Lacs analysed in CG08 
and with the FSRQs analysed in Maraschi et al. (2008, hereafter 
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Log [Hz] Log [Hz] 



Figure 4. SED of 0917+449=84 0920+44, 1013-054 and 1329-049. Sym- 
bols and lines as in Fig.[T] For PKS 0917+449 we show two models, corre- 
sponding to the Fermi/Swift data and to the older EGRET data. 



Figure 5. SED of 1424-354, PKS 1502+106; PKS 1520+319 Symbols and 
lines as in Fig.[T] 



MOS). As can be seen, the blazars analysed in CG08 and M08 tend 
to have larger jet powers with respect to disk luminosities than the 
blazars presented here. This can be partly understood considering 
that a large fraction of the blazars in CG08 were EGRET sources, 
and the 7-ray data correspond to high state of the sources. In our 
sample, instead, the 7-ray luminosity, derived by the average 7- 
ray flux seen by Fermi, is considerably smaller than the maximum 
values seen by EGRET (see Fig. [8] as illustrative example). The 
smaller electromagnetic power output implies that also the kinetic 
power carried by the jet is smaller The blazars considered in 1VI08 
(see also Sambruna et al. 2007 and Ghisellini 2009), instead, lack 



7-ray data, and the estimates were based on fitting the SED up 
to the Swift-XKI energy range. Several of them showed a very 
hard X-ray spectrum, suggesting a peak in the 1-10 MeV region 
of the spectrum. IVIeV blazars of this kind, peaking at energies sig- 
nificantly smaller than 100 MeV, may not be easily observable by 
Fermi, and could be among the most luminous blazarfl (according 
to the blazar sequence, Fossati et al. 1998). The jet powers calcu- 
lated here may not correspond to the "very high states" of blazars as 

^ The recently published list of blazars detected in the 15-55 keV band in 
the 3-years survey of the Burst Alert Telescope (BAT) onboard Swift indeed 
includes very powerful blazars (Ajello et al. 2009), up to redshifts z = A.l. 
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Figure 6. SED of PKS 1551+130, PKS 1633+382 and PKS 2023-077. 
Symbols and lines as in Fig.[T] 
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Figure 7. SED of PKS 2052^7, PKS 2227-088, PKS 2052^47 and PKS 
2325+093. Symbols and lines as in Fig.[T] 



seen by EGRET. They are therefore more representative of a typical 
blazar state, but we have to remember that when constructing a list 
of variable sources (even with a more sensitive instrument) high 
states are always over-represented. What is confirmed with these 
new data is that, at least occasionally, the jet is more powerful than 
id. 

We can also compare the fraction of the total jet power carried 
in radiation (er), magnetic field (eb), and electrons (tc) with the 
sample of blazars studied in CG08. This is done in Fig.[T3] showing 
a rough consistency with the blazars in CG08. 



6.3 Emission mechanisms 

The location of the dissipation region is almost always within the 
BLR (exceptions are 0215+015 and 1520+319), that provides most 
of the seed photons scattered at high frequencies. The main emis- 
sion process are then synchrotron and thermal emission from the 
accretion disk for the low frequency parts, and External Compton 
for the hard X-rays and the 7-ray part of the spectrum. The X-ray 
corona and the SSC flux marginally contribute to the soft X-ray 
part of the SED. The IR radiation reprocessed by the torus can in 
some case dominate the far IR flux, but the scarcity of data in this 
frequency window does not allow to fully assess its importance. 
Note that when i?diss < ^^blr, the overall non-thermal emission 
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Figure 10. Histograms of the total jet power and the accretion disk lumi- 
nosity, both in cgs units (left) and in Eddington units (light). 
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Figure 11. Histograms of the outflowing (top left) mass rate; the accretion 
mass rate (bottom left) and their ratio (top right). On the bottom right we 
show the distribution of black hole masses found in this paper. 
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Table 5. Jet power in the form of radiation, Poynting flux, bulk motion 
of electrons and protons (assuming one proton per emitting electron). Sub- 
scripts refer to the year of the corresponding SED or to an EGRET (_E) state 
as opposed to a "Fermi" (F) state. 



is rather insensitive to the presence/absence of the IR torus, since 
the bulk of the seed photons are provided by the broad lines. How- 
ever, for 0215-1-015 and 1520-1-319, the location of the dissipation 
region is beyond the BLR, and for these two sources the exter- 
nal Compton process with the IR radiation of the torus is crucial 
to explain their SED. Since their dissipation region is beyond the 
BLR, the corresponding energy densities in radiation and fields are 
smaller than in the other sources, as can be seen in the top panel of 
Fig. [14] Furthermore, the smaller Uy + Ub implies a slower cooling, 
and thus a larger 7cooi (bottom panel of Fig.ll4b. 



6.4 Particle distribution 

The lower end of the emitting particle distribution is usually the 
bottleneck to reliably calculate how many particle the jet carries, 
and hence its power. This is because low energy electrons emit, 
by synchrotron, in the self-absorbed, not visible, regime, and, by 
SSC, they contribute in a region of the SED usually dominated by 
synchrotron or by the accretion disk. There are however several 
cases where the X-ray data can be convincingly fitted by the EC 
mechanism, and in these cases the IC radiation produced by these 
low energy electrons is dominating the observed flux. These are 
the cases where the X-ray spectral shape is very hard, especially 
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at lower frequencies, excluding SSC as the main contributor to the 
observed flux. As examples, see the BeppoSAX spectrum of 225 1- 
158 in the year 2000, or the XRT spectrum of 2325+093. These 
are the cases where the A''(7) distribution can be reliably deter- 
mined even in its low energy part. Theoretically, as an outcome of 
the model assumption, we can also calculate the random Lorentz 
factor of the electrons cooling in one light crossing time. This de- 
pends on the amount of radiative cooling, in turn depending on the 
bulk Lorentz factor, the magnetic field, and the amount of exter- 
nal radiation. Fig. [14] shows the values of 7cooi (as a function of 
the peak energy 7pcak) for our sources. Since we are dealing with 
very powerful blazars, where the cooling is strong, the derived val- 
ues are particularly small, of the order of 1-10. This ensure that 
Nl'y) (X between 7cooi and 7b. 

For all our sources the importance of the 7-7 process 
(which is included in our model) is very modest, and does not in- 
fluence the observed spectrum. 



6.5 Comparison with low power jets 

Fig. [m shows the value of the peak energy 7pcak (namely the ran- 
dom Lorentz factor of the electrons producing most of the emission 
at the two peaks of the SED) as a function of the total energy den- 
sity (as seen in the comoving frame; upper panel) and as a function 
of 7cooi (calculated after one light crossing time rdiss/c; bottom 
panel). We have reported, for comparisons, the values obtained by 
CG08 for a sample of blazars including both powerful FSRQs and 
low power BL Lacs. As expected, the values for our blazar are con- 
sistent with those previously found by CG08 for powerful FSRQs. 
The bottom panel of Fig.[T4]is a different version of the blazar se- 
quence as originally proposed in Ghisellini et al. (1998) and shown 
in the upper panel. It suggests that 7pcak becomes equal to 7cooi 
when the cooling is weak, but it must be larger than 7cooi when 



Figure 13. The fraction of Ljct radiated (tr, top panel), in relativistic lep- 
tons (to, mid panel) and in magnetic fields {eg, bottom panel) as functions 
of Pjct = Pp + Pe + Pb- Our sources (black circles) are compared to 
the FSRQs (grey triangles), BL Lacs (grey empty circles) and TeV BL Lacs 
(grey filled circles) of the blazar sample studied in CG8. 



the cooling is fast. Therefore the position of the peaks of the SED 
of strong FSRQs is controlled by the injection function Q{j) and 
its break at 7b, that is proportional to 7pcak (the two energies are 
not exactly equal because of the assumed smoothly curved injec- 
tion function, see GT09). In the slow cooling regime, instead (low 
power BL Lacs), 7b becomes unimportant, and the particle distri- 
bution self-adjust itself. The new feature added here is the fact that 
in this case (i.e. slow cooling) the value of 7pcak ~ 7cooi depends 
also on the size of the source, and not only on the total energy den- 
sity. 

6.6 Multi-states of specific sources 

PKS 0528+134 — Fig. [3] shows 4 states of this source, one of the 
best studied by EGRET, that was able to detect it at very different 
levels of 7-ray activity (see also Mukherjee et al. 1999 for a de- 
tailed analysis of the different states during the EGRET era). Keep- 
ing the black hole mass and the accretion rate fixed (so that also 
Ld and -Rblh are fixed), we tried to interpret the 4 different states 
changing a minimum number of parameters. Also the bulk Lorentz 
factor need not to change. Indeed, the main changing parameters 
are i) i?diss, implying a change of the magnetic field in the dissipa- 
tion region (larger field at closer distances, according to a constant 
Pb, as listed in Tab.[5]l, and ii) the power P( injected into the source 
in the form of relativistic electrons. The change in P( is the main 
responsible for the changes in the jet bulk kinetic power (factor 
^5), that are sufficient to explain the rather dramatic changes in 
the observed SED. We had to change also the high energy slope 
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sources with i?diss > ^?BLR and this impUes a reduced energy density (in 
radiation and magnetic field) and therefore a larger 7cool- 



(S2) of the injected electron distribution, and (but by a modest fac- 
tor) 7i and 7max- The changes of the magnetic field (factor ~2) 
are responsible for a very different contribution of the SSC flux in 
the X-ray band, dominated by this component in the SEDs of 1993 
and 1995. 

TXS 0917+449 = RGB 0920+446 — In Fig. g] we show 2 models 
for the SED of this source, corresponding to the combined Swift 
and Fermi data, and also to the old EGRET ones (for which there 
is a paucity of data). Therefore the two models do not describe two 
different observed SEDs, but can illustrate how the model param- 
eters have to change if only the 7-ray emission changes. In this 
case we have fixed, besides the black hole mass and accretion rate, 
also the location of the dissipation region, the magnetic field and 
the injected power We only changed the high energy slope of the 
injected electron distribution (S2, from 2.6 to 2), and 71 (from 50 
to 30). This changes are responsible for the (small) changes of the 
total bulk kinetic power of the jet. 

PKS 1454-354 — Fig.|5]shows 3 SEDs for this source, correspond- 
ing to three states already observed by Fermi and discussed in Abdo 
et al. (2009c). Fig.[5]reports the 7-ray data of the flare, of the "low" 
state and of the intermediate 7-ray state that corresponds to the av- 
erage flux reported in Abdo et al. (2009a). Swift AsXa corresponds to 
Jan and Sept 2008. As for the other sources, we have tried to model 
the 3 states of the source assuming the same black hole mass and 
disk luminosity (hence the same -Rblr)- Values of i?diss are con- 
strained by the observed short variability timescale of the 7-ray 



flux (Abdo et al. 2009c). We can reproduce the different source 
states by changing JJdiss (factor 2), Pi (factor 5), F (almost a fac- 
tor 2) and B (factor ~2.5). In addition, we had to adjust 71 and S2. 
These changes of the physical parameters of the jet imply differ- 
ent jet powers (factor ~2 for Pb and ~10 for the bulk kinetic jet 
power). 

2251-158 = 3C 454.3 — Fig.|9]shows 4 states of this source, one of 
the brightest AGILE blazars, and also detected, besides by Fermi, 
also by EGRET. The first 3 SEDs have been discussed by Ghis- 
ellini et al. (2007), with the difference that at that time only an inte- 
grated AGILE flux was known, while now (Vercellone et al. 2009) 
the spectral analysis of the AGILE data yielded a spectrum rising in 
I'Fi,. The fourth corresponds to the Fermi observation (with a flux 
slightly smaller than AGILE but with a steep slope). This blazars 
shows a rather dramatic variability at all frequencies, most notably 
in the optical band in 1995 (Fuhrmann et al. 2006; Villata et al. 
2006; Giommi et al. 2006; Plan et al. 2006). Comparing the SED 
of the 2005 flare with the SED of 2000, the IR-optical flux in- 
creased by two orders of magnitude, while the hard X-rays flux 
increased by a a factor ~10. Unfortunately, for the SED in 2005 
there are no information about the 7-ray flux. Katarzynski & Ghis- 
ellini (2007) interpreted this behaviour (namely, more variability 
amplitude in the optical than at X-rays) suggesting that the main 
change was due to the location of the dissipation region in the jet, 
more internal in 2005. They showed that even with a constant jet 
power, and a constant bolometric luminosity, dramatic variations 
in specific bands were possible. The reason is that a more internal 
dissipation region may correspond to a larger magnetic field and a 
smaller bulk Lorentz factor (i.e. the jet may still be accelerating): 
the combination of these two key factors yields a larger magnetic 
energy density and a smaller external radiation field (since in the 
comoving frame it is less amplified by the smaller F). This in turn 
implies more synchrotron and less EC radiation. These ideas were 
later confirmed (Ghisellini et al. 2007) by the first AGILE obser- 
vations, that caught the source in a bright 7-ray state and in a rel- 
atively faint optical state (compared to 1995). We here follow the 
same ideas, adopting a smaller i?diss and F for the 2005 state (see 
Tab. |4j. The SED corresponding to the Fermi/Swift observations 
in 2008 confirms this scenario: the source is still much brighter in 
7-rays than revealed by EGRET, but it is not exceptionally strong 
in optical. What is notably different, for the 2007 and 2008 states, 
is the slope of the 7-ray spectrum, which is flat (in vF^) for AG- 
ILE and steep for Fermi. We modelled the different 2007 and 2008 
states with a different injected electron distribution, much flatter in 
2007, leaving all other parameters unchanged. This implied a dif- 
ferent number of emitting electrons, and a somewhat different jet 
power (see Tab.|5). 



7 DISCUSSION 

Several of our blazars have good Swift data that allowed to estimate 
the black hole mass and the disk luminosity. There are two reasons 
for the disk luminosity to emerge in the optical UV range of the 
spectrum in these luminous sources: i) according to the blazar se- 
quence, and directly seen in the Fermi data, the high energy spec- 
trum is steep, implying a high energy peak below the Fermi en- 
ergy range. This implies that the corresponding synchrotron spec- 
trum peaks in the IR band, with a relatively smaller contribution 
in the optical-UV. ii) If highly luminous blazars have large black 
hole masses, then we expect relatively smaller magnetic fields in 
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the dissipation region with respect to blazars with smaller masses 
(see Ghisellini & Tavecchio 2009), and thus an enhanced external 
Compton component with respect to the synchrotron one. Since 
ours are the most luminous Fermi blazars, is likely that they have 
black holes with the very large masses. This is is indeed what we 
found: all (but one) of them are above 10^ Mq. We can compare 
these values with the ones found, with a different method, for the 
quasars found by the Sloan Digital Sky Survey (SDSS) (see e.g. 
Vestergaard & Osmer 2009; Vestergaard et al. 2008; Kelly, Vester- 
gaard & Fan 2008; Natarajan & Treister 2009; Shen et al. 2008). 
The method to find the black hole mass in this case is to measure 
the width of some broad emission line and the luminosity of the op- 
ticalAJV continuum. Since the latter is found to correlate with the 
typical size of the BLR (e.g. Kaspi et al. 2007; Bentz et al. 2006), 
the hypothesis of virial motion of the broad line clouds allows to 
estimate the black hole mass. The corresponding uncertainties are 
large, but studying very large samples of thousands of quasars it 
emerged that there is a "saturation" value of the black hole mass, 
around 5 x lO^-lO^'Af©. Since our FSRQs are among the most lu- 
minous blazars, their black hole masses are in agreement with what 
found by the SDSS. This consistency, found using two different 
methods for the mass determination, is reassuring. 

We expect that a powerful jet is accompanied by a close-to- 
Eddington accretion rate. This is based on the relation between the 
jet power and the accretion disk luminosity, as pointed out in sev- 
eral works. The two powers are found to be approximately equal in 
Rawlings & Saunders (1991). More recent work showed that jets 
are ~10 times more powerful than the luminosity of their disks 
(CG08; Maraschi & Tavecchio 2003). 

By model fitting, we could derive the intrinsic physical quan- 
tities needed to estimate the power carried by the jet. Therefore we 
confirm the previous hints of a link between the jet power and the 
disk luminosity. On average, the jet power is larger (but it is of the 
same order of) the disk luminosity, suggesting Pj ~ Mc^ . For sev- 
eral of our blazars the VWOT/Swift data allowed a direct estimate of 
the disk luminosity, that together with the black hole mass estimate 
give us the disk luminosities in Eddington units: the values clus- 
ter around L^/L-Edd ~ 0.1 — 0.3. We re-iterate that these results 
are very approximate and must be taken with caution, given the 
uncertainties involved (especially in decomposing the non-thermal 
and thermal fluxes, and the possibility to have some residual dust- 
extinction in the host fame). One should also take into account the 
possibility that the observed disk luminosity never reaches the Ed- 
dington value in standard disks, since some of the produced radi- 
ation may be trapped within the accretion flow and swallowed by 
the black hole before it can be radiated away. 

Our sample, by construction, includes only very powerful 7- 
ray blazars: for these sources it is more likely that the accretion 
disk directly shows up in the SED. However, this disk component 
is seen also in blazar with not so extreme luminosities as in 3C 
273, where it is always visible, and in PKS 1510-089 (Kataoka et 
al. 2008) and 3C 279 (Pian et al. 1999). In Fig.[T2]we have marked 
the positions of these sources in the Pj-Ld plane. There might be 
different reasons for the different "visibility" of the accretion disk 
in blazars. In 3C 273 its luminosity is of the order of 3 x 10^" er 
s~^, of the same order than the non-thermal observed luminosity. 
In Celotti & Ghisellini (2008) we found that the Doppler boosting 
in this source was slightly less than for the rest of the blazars we 
analysed. We therefore suggest that for 3C 273 (a relatively nearby 
source, z = 0.158) we are looking at a jet very slightly misaligned 
(say 5° instead of 3°), letting the disk emission to always stick out 
from the non-thermal spectrum. 



For the other two sources the accretion disk luminosity is 
much less than in 3C 273 (i.e. of the order of 10*^ erg s~^). This 
can be the consequence of a combination of a smaller black hole 
mass and a smaller accretion rate with respect to 3C 273, leading 
also to a somewhat smaller jet power (see Fig.ll2t. 

Since almost all dissipation regions are located within the 
BLR, the most relevant seed photons for scattering are Lya pho- 
tons. As explained in Ghisellini & Tavecchio (2009), this implies 
a break in the intrinsic spectrum at~ 15/(1 + 2:) GeV due to the 
decreasing (with energy) Klein-Nishina cross section and to the 
7-7 process between high energy and broad line photons. 

If spectra will not cut-off at 15/(1 + z) GeV, the seed photons 
must be at lower frequencies (Ghisellini & Tavecchio 2009), as oc- 
curs if the dissipation region is beyond the BLR, where the most 
relevant seed photons becomes the IR photons reprocessed by the 
torus. In this case we also expect larger sizes of the emitting blob 
(since it is located at larger distances from the black hole), hence a 
relatively slower variability. As pointed out in Ghisellini & Tavec- 
chio (2009), these are also best candidates to study the effect of the 
photon-photon absorption process due to the intergalactic optical- 
UV backgrounds. 
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